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We study large-amplitude one- dimensional solitary waves in photonic crystals featuring compe- 
tition between linear and nonlinear lattices, with minima of the linear potential coinciding with 
', maxima of the nonlinear pseudopotential, and vice versa {inverted nonlinear photonic crystals, IN- 

(S| ■ PhCs), in the case of the saturable self- focusing nonlinearity. Such crystals were recently fabricated 

using a mixture of SU-8 and Rhodamine-B optical materials. By means of numerical methods and 
analytical approximations, we find that large-amplitude solitons are broad sharply localized stable 
pulses (quasi-compactons, QCs). With the increase of the total power, P, the QC's centroid performs 
multiple switchings between minima and maxima of the linear potential. Unlike cubic INPhCs, the 
^3 ' large-amplitude solitons are mobile in the medium with the saturable nonlinearity. The threshold 

CO , value of the kick necessary to set the soliton in motion is found as a function of P. Collisions 

between moving QCs are considered too 

Q ; PACS numbers: 42.65.Tg; 42.70.Qs;05.45.Yv 
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I. INTRODUCTION AND THE MODEL 



Photonic and matter waves propagating under the combined action of linear and nonlinear lattice (LL and NL) 
\ potentials exhibit a plenty of dynamics In particular, solitons in such media is a topic of considerable current 
interest, see recent review [2^ and original works [1,13 devoted to the studies of one-dimensional (ID) solitons. Other 
works were dealing with solitons [H, [gI and localized vortices 0] in two-dimensional(2D) versions of such systems, 
which represent, in particular, photonic crystal fibers [H and 21) photonic crystals [9| made of nonlinear materials. 
^ ■ The evolution of solitons in these systems obeys the nonlinear Schrodinger equation (NLSE), in which the LL and NL 
On : are represented, respectively, by a usual periodic potential and by a periodic pseudopotential [H, [IqI, which is induced 
^ ■ by a spatially periodic modulation of the local nonlinearity coefficient. In optics, the NL represents the mismatch 
between the nonlinearity of the host material and the stuff filling voids of the photonic-crystal-fibers structure, which 
may be air, another solid material or a liquid crystal p^]. The same model, in the form of the Gross-Pitaevskii 
equation, applies to matter waves in a Bose-Einstein condensate (BEG) which is trapped in a combination of a linear 
periodic potential, created by an optical [l3| or magnetic lattice, and a pseudopotential lattice, that may be 
induced by a periodic modulation of the local nonlinearity provided by a properly patterned [2| external magnetic 
[m or optical fields [Ts]. 

Recently, a combination of competing 7r-out-of-phase-juxtaposed LL and NL, with maxima of the refractive index 
coinci ding with minima of the local strength of the self-focusing nonlinearity and vice versa, was considered in several 
^ ■ works pj|-[l9|. This medium may be naturally called an inverted nonlinear photonic crystal (INPhC). It was reported 
I that INPhCs could be fabricated by means of the technique based on direct laser writing in silica [20]. It was also 
predicted that a similar setting can be created in a virtual form, using the electromagnetically induced transparency 
acting on dopant atoms periodically distributed in a passive matrix [l9|. Due to the competition between the LL 
and NL, solitons in INPhCs may feature specific power-dependent properties, such as double symmetry breaking [l9|. 
While previous studies of INPhCs were dealing with the Kerr (cubic) nonlinearity, in this work we consider solitary 
waves in the system with saturable nonlinearity, which occurs in various optical media. The extension of the analysis 
for this nonlinearity is natural, whilst solitons in INPhCs feature the strong sensitivity to the power. 

First, we propose an experiment setup to realize an INPhC with a saturable nonlinearity. Recent works reported 
the creation of resonantly absorbing waveguide arrays and imaginary-part photonic crystals, that feature a spatially 
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FIG. 1: (Color online) (a) The one- dimensional nonlinear photonic crystal, with the blue and gray areas depicting the nonlinear 
and linear stripes, respectively, (b) The wavelength dependence of the absorption coefficient and refractive- index variation, a 
and An (solid and dashed lines, respectively) in the SU-8 - RhB mixture. 



periodic modulation of the absorption coefficient, built on the basis of the SU-8 polymer (a commonly used transparent 
negative photoresist) doped with Rhodamine B (RhB, a dye featuring saturable absorption) [111, [22|, see Fig. [T]for a 
schematic setup. 

A similar pattern can be used for our purposes. In the paraxial approximation, the light propagation in the array 
obeys the spatial-domain NLSE for the local amplitude of the electromagnetic field, u{z^x)^ where z and x are the 
propagation distance and transverse coordinate: 

1 {kQ8n^ia/2)R{x) 

lUz = —ITtUxx (1) 

Here, R{x) is the array's structure function which is to be taken as per Fig. |l(a)[ k — k^n^ where n = 1.62 
is the refractive index of pure SU-8, ko = 27r/A, and A is the wavelength of light in vacuum. Further, a and 
Sn are the absorption coefficient and refractive index difference between the SU-8 - RhB mixture and pure SU- 



8, respectively, which are coupled by the Kramers-Kronig relation, and are plotted versus A in Fig. 1(b) [22|. 
According to experimental data [22), in the relevant range of far blue shift from the absorption peak (A ^ 460 nm), 
a/2 (l/20)/co \Sn\^ which allows one to neglect the dissipative term in Eq. ([1]), simplifying it to 



1 y{x) 

lUz = — Uxx H r-i^u, (2) 

2 1 + \u\^ ' ^ ^ 

where k was removed by rescaling of x, and V{x) = — {koSn) R{x). 

In the system under consideration, light experiences strong saturable self-focusing induced by RhB in the mixed 
material, while the nonlinearity of pure SU-8 is negligible. On the other hand, the refractive index of the mixture is 
smaller than in pure SU-8, i.e., (5n < at A = 460 nm, hence function V{x) in Eq. (|2]) takes positive values, and the 
waveguide array meets the condition of the tt phase shift between the spatial modulations on the linear and nonlinear 
local characteristics of the medium, thus realizing the INPhC with the saturable nonlinearity. 

The objective of this work is to study the existence, stability, mobility, and interactions of solitons in the INPhC 
model based on Eq. (j2j). In Sec. II, we study properties of the solitons by means of numerical simulations and 
analytical approximations. It will be demonstrated that they feature sharp localization, i.e., a quasi-compacton (QC) 
shape. Following the increase of the total power, P, the QC switches its position between maxima and minima of the 
linear refractive index. In Sec. Ill, we study mobility of the QCs, imposing the phase tilt onto them, i.e., suddenly 
multiplying the wave form by exp(i7^x). A critical tilt (alias kick)^ Vc{P)^ beyond which the compacton starts to move, 
is found. The dependence r]c{P) features variations with the same period in P as the above-mentioned switching of the 
quiescent solitons. Collisions between moving compactons are also studied in Sec. Ill by means of direct simulations. 
The paper is concluded by Sec. IV. 
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(d) (e) (f) 

FIG. 2: (Color online) Numerically found soliton profiles (blue solid lines) for the modulation function (shown by the black 
solid line) with depth Vb = 0.02 (recall the modulation period is scaled to be d = 20). (a) At total power P — 400, the center of 
the soliton is located at the minimum of V{x). (b) At P — 800, the soliton shifts to the maximum of V{x). (c) At P — 2000, the 
solitons shifts back to the minimum of V{x). Panels (d), (e) and (f) illustrate the stability of the solitons via direct simulations 
of their perturbed evolution. 



II. QUASI-COMPACTONS: NUMERICAL AND ANALYTICAL RESULTS 



A. Numerical simulations 



Modulation function V{x) in Eq. (|2]) corresponding to Fig. [T]is a piecewise-constant one, of the Kronig-Penny type. 
In this paper, we approximate V{x) by the first term of its harmonic decomposition, assuming that contributions of 
higher harmonics are negligible for sufficiently broad solitons: 

V{x) = (Fo/2) [1 - cos {271x1 d)] , (3) 

where d and Vb > are the modulation period and depth. The scaling is set by fixing d = 20, which leaves Vq as a 
free parameter. In this section, we present numerical results for Vq = 0.02, 0.03 and 0.04, which adequately illustrate 
the generic situation. Below, we focus on solitons with a sufficiently large amplitude, as for small amplitudes the 
truncated expansion of the saturable nonlinearity amounts to the previously studied Kerr model pjl-IToj. 

Typical numerical results for solitons solution at different values of the total power, P = \u{x)\'^dx^ are 

displayed in Fig. [2l The stationary profiles, presented in Figs. 2(a)|2(c)^ were generated by dint of the the imaginary- 



time-propagation method [23|. It is concluded that the profiles are strongly localized, making the solitons QCs 
(quasi-compactons) [a solution fully localized within a finite (compact) interval of y is usually called a compacton 
05 mi]- On the other hand, the soliton broadens with the increase of P. At relatively small values of the total power, 
e.g., P = 400, the center of the soliton is located at a minimum of V{x) [Fig. [2(a)] ; then, at P = 800 [Fig. [2(b)] , the 
position of the soliton is switched to a maximum of V{x)^ and at P = 2000 it switches back to the minimum [Fig. 

To quantify power-dependent properties of the solitons, we define the center-of-mass coordinate and average width 
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FIG. 3: (Color online) The plots in panels (al)-(cl) display the coordinate of the soliton's center of mass, Xmc, as a function 
of P, for Vo = 0.02, 0.03 and 0.04, respectively. (a2)-(c2). The critical kick which initiates the motion of the soliton, ?7c, as a 
function of P for Vo = 0.02, 0.03 and 0.04. (d)-(f) The average width of the soliton as a function of P for Vo = 0.02, 0.03 and 
0.04. 



of the soliton as follows: 



/ + 00 
x\U{x,P)\''dx 
-oo 

/ + 00 
{x-X^ef\U{x,P)\'- 
-oo 



(4) 



where U (x, P) is the soliton solution with power P. Using the integration in imaginary time, the data were collected 
for the solitons from P = 100 to 10000, with modulation depths Vq = 0.02, 0.03 and 0.04. Real-time simulations of 
the evolution of these solutions in the framework of Eq. (j2j) demonstrate that they all are stable against perturbations. 

Numerically found dependences dH are displayed in Fig. [3l The upper panels. Fig. [3l^a-c), show that the soliton's 
center of mass switches (as said above) between the minimum {x = 0) and maximum (x = —10) of the modulation 
function V{x). The period of the switching, AP, increases with the increase of modulation depth Vq. Further, Figs. 
3(d)||2(f)| show that the average width of the soliton grows with the increase of P and decrease of Vq. 



B. Analytical approximations 

In saturable nonlinear media, large-amplitude solitons are always broad. Accordingly, to develop an analytical 
approximation, we assume a soliton with the width much larger than d. Then, in the zero-order approximation, Eq. 
(|2j) is replaced by the equation with the averaged potential, V{x) 

Equation (|5]) gives rise to soliton solutions, u{z,x) = e^^^U{x — ^), with the coordinate of the soliton's center ^, and 
real function U (y) obeying the following equation: 
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where y = x — The effective potential in Eq. (j6j), 

mf/) = fc/7^ + |ln([^), (7) 

formany corresponds to the Newton's equation of motion for a particle with coordinate U (y) {y plays the role of time) 
and unitary mass. Then, the shape of the soliton is determined by the respective energy equation, 

(f)%2.W-C/') + ,,,„(i±M)=„, (8, 

where Uq is the amplitude of the soliton (the largest value of U). As usual, the soliton trajectory corresponds to the 
solution of Eq. (jS]) starting from U = Oat?/ = — oo, bouncing back from the potential well at[/ = [/o, ^ = 0, and 
returning to U = Oat^^+oo. Setting in Eq. (j8j) = dU/dy = 0, one can find a relation between k and Uq. Being 
interested in solutions with large amplitudes, we assume here /7o ^ 1, which yields 

k^-{Vo/2)Uo'ln{Ui), (9) 



f^min«t^o/ln(C/2), (10) 

Umin being the coordinate of the minimum of potential ([7j). Equation (p!Q|) implies [/^-^ <C ^/q , i.e., the minimum of 
the potential is located much closer to /7 = than the largest value Uq. This fact suggest a possibility to use the 
following approximation for solving Eq. ([8j): at the first stage of the approximation, we drop the logarithmic term in 
Eq. ([8|) altogether, i.e., we replace the equation by its "primitive version", 

^)\2k{Ui-U^)=0, (11) 
whose obvious solution is 



U{y) = Uo cos (v^y) « Uo cos |^ V%MMy j , (12) 

where Eq. (|9]) was used to replace k by the Uq. 

At the second stage of the approximation, we recall that the particle does not perform periodic oscillations, as 
formally follows from the expression (p!2]) . but it starts the motion from U = Oat?/ = — oo, and returns to U = at 
y +00, as the soliton solution must do, see above. This means that approximation (fT2]) is usable in interval 

When the moving particle approaches edges of interval (p!3|) . the logarithmic term in Eq. (j8j), which was neglected in 
Eq. (fTT]) . "suddenly" becomes important, leading to the stoppage of the particle. Thus, the full approximation in the 
large-amplitude limit, /7o ^ 1, amounts to representing the soliton as the QC (quasi-compacton): 

Uiy) ^ { ^^""""^ I^ A^M^o) ^^ ^ ^ ^^^^^2, ^^^^ 

0, at 1^1 > Ltot/2. 

where Ltot, given by Eq. ([13]), plays the role of the total width of the QC. Note that the large amplitude Uq of the 
soliton implies that width ([13]) is also large, i.e., the QC is a sharply locahzed but wide localized pattern. In the 
present approximation, the total power of the soliton can be calculated as follows: 



+00 1 tt3 



Further, using the fact that P and Uq are large, Eq. ([15]) can be approximately inverted, to give the peak power 
(squared amplitude) of the QC as a function of its total power: 

— \ 2/3 
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(a) (b) 

FIG. 4: (Color online) (a) The comparison of the numerical results for the peak power of the quasi-compactons with the 
analytical prediction given by Eq. ()16p . (b) The full width of the quasi-compacton, as predicted by Eq. ()17p . In this figure, 
Vq = 0.02. 



Next, using Eq. (fT3|) . it is also possible to approximately express the full width of the compacton in terms of the total 
power: 



3^" P 
Vq InP 



1/3 



(17) 



Relations (p!6]) and ([TTj) are plotted in Fig. |4]for different values of P. In particular, Fig. |4l^a) demonstrates that 
the analytical prediction for the peak power is consistent with the numerical results. On the other hand, the values 
of the width given by Eq. (pT|) [see Fig. Hl^b)] are larger than their numerically found counterparts presented in Fig. 
[3jd). However, one should take into regard that width ([17]) is "all-inclusive" (total), while the width shown in[3l^b) 
is the average one (definitely far smaller than the total width). Note also the similarity in the P-dependence of the 
width in Figs. [S^d) andg^b). 

To explain the position switchings of the QC, one may use the effective energy (potential) of the interaction of 
the QC with the spatial modulation, corresponding to the term ~ cos(27rx/(i) in expression ([3]). This energy can be 
defined as follows, with regard to the definition of ^ = x — ^ in Eq. (|6]): 



EUO = -Vo J dxcos (^^) In [l + U\x)] 

pi+Ltot/2 /27rr\ 
-Vo / cos — In [1 + U^{x)] . 



/{-Ltot/2 



(18) 



In this expression, it is taken into account that QC (p!4|) occupies, approximately, a finite region of x, ^ — Ptot/2 < 
X < ^ + Ltot/2. After three integrations by parts, expression ([18]) can be approximately calculated as follows. 



cos 



1^ 



where = V^d^/n^. Expression ([T9)) predicts that the minimum of the interaction energy is located at points ^ 
with integer n (at minima of the modulation function) for values of Uq such that sin (^ir'^Uo/ 

and at points = (n + 1/2) d (at maxima of the modulation) in the opposite case, sin (^tt'^Uq/ 



dv/K,ln(C/2) 



(19) 



= nd 
< 0, 



> 0. 



Because Uq grows as the function of P, as per Eq. ()T6|) . this means that, with the increase of P, the position of the 
energy minimum must indeed switch between adjacent extrema of the modulation function. The fact that Vq appears 
in the denominator of the argument of the sinusoidal factor explains why the period of the switching increases with 
Vo. 



7 





(a) 



(b) 



FIG. 5: (Color online) Direct simulations of the compacton with P = 2000, initiated by the application of the phase tilt with 
r] = 0.020(7r/d) (a) and r] = 0.025(7r/(i) (b). The lattice depth is Vo = 0.02. 



III. MOBILITY OF THE QUASI-COMPACTONS 
A. Driving the soliton by the phase tilt 



Mobility of solitons under the combined action of LL and NL is a problem which was considered in a number of 
different settings pi, l26l - [29| . To set the QC in motion, we follow the standard approach, suddenly kicking a quiescent 
one, i.e., multiplying the respective solution, /7(x), by the phase-tilt factor: U{x^r]) ex.-p{ir]x)U{x)^ where r] is the 
strength of kick. 

The simulated propagation of the QC, induced by kicking the standing one, Uix^P = 2000), is displayed in Fig. [51 
for two different values of rj. It is observed that r] = 0.020(7r/(i) causes only oscillations of the soliton, without any 
progressive motion. However, a slightly higher kick, r] = 0.025 (yr/d), is sufficient to set the same soliton in persistent 
motion. These results imply that there must be a threshold (critical) value of the kick, 77c, such that the kick with 
T] > r]c makes the given QC a traveling soliton. For P = 2000, it is rjc = 0.021 {Tv/d). 

The critical kick for the QCs is shown as a function of the total power in Figs. [3l^a2)-|3l^c2), for Vq = 0.02, 0.03, and 
0.04. The general trend to the decrease of rjc with the increase of P is explained by the broadening of the soliton with 
the increase of the total power in the case of the saturable nonlinear ity, and ensuing enhancement of the solitons' 
mobility [29^. A specific feature of the INPhC is that oscillations of r]c are superimposed on top of the gradual decay 
of the critical kick, the period of the oscillations coinciding with that of the switching of the center-of-mass position. 
Figures 3(a)||3(c) also demonstrate that rjc d ecreases w ith the decreases of modulation depth Vq. The latter trend is 



correlated with the fact that, as seen in Figs. [3(d)p(f)| the width of the solitons increases, making them more mobile, 
with the decrease of Vq (as well as with the increase of P). 

Further, Fig. [6] shows that the momentum of the moving soliton, which is defined as 

{p{z)) = -iJ dxU{x,z,r]yU,{x,z,r]), (20) 

is, naturally, a linear function of the initial kick, 77. Further simulations (not shown here) demonstrate that the QC 
are not destroyed even by a very strong kick, far exceeding rjc 

It is relevant to compare the mobility of large- amplitude solitons in the INPhCs models with the saturable and 
cubic self-focusing. In the latter case, Eq. (|2]) is replaced by Q 

iuz = -^u^x + V{x){l - \u\^)u, (21) 

where we chose the modulation function, V{x)^ in the same form as in Eq. ([3]), with Vq = 0.02 and d = 20. A 
typical example of simulations of the mobility of solitons in this model is displayed in Figs. [71^a)-(c), for three kicks, 
T] = 0.157r/(i, 0.257r/(i, and O.Sbn/d. In particular. Fig. |7(a)| sho ws that rj = O.IStt/^ give rise to oscillations of the 
soliton, without depinning. Further, it is observed in Fig. 7(b) [ that the stronger kick, r] = 0.257r /(i, tr ansforms the 



soliton into an oscillatory beam with a fuzzy shape, but, still, it was not set in motion. Finally, Fig. 7(c) [ demonstrates 



that the strongest kick, 77 = O.Sbn/d^ completely destroys the soliton (recall the destruction never happens in the 
model with the saturable nonlinearity ) . Therefore, there must exist critical values of the kick separating the stable 
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FIG. 6: The momentum of the sohton with P = 2000, measured dut z = 1000, as a function of initial kick 77. 



solitons, fuzzy beams, and complete destruction. These borders are plotted in the top panel of Fig. |7(d)[ in the range 
of powers 50 < P < 100 (the average width of the static soliton in the same region is plotted, versus P, in the bottom 
panel). Thus, the large-amplitude solitons in the INPhC model with the cubic self- focusing are not mobile at all. On 
the other hand, for dynamics of small-amplitude solitons in the models with the saturable and cubic nonlinearities 
are similar, due to the obvious expansion, (l + \u\'^) ^ 1 — \u\'^ (results for this case are not shown here in detail, 
as they are less interesting). 



B. Collisions between moving quasi-compactons 



Collisions is a natural way to test interactions between solitons. Here we present some numerical results for collisions 
between QCs in the INPhC with the saturable nonlinearity. The initial state at 2: = is taken as a set of two far 
separated kicked solitons: 



u{x) = /7(x + xo,P)e^^^ 



U{x — xo, P)e" 



-ir]2{x — xo) 



(22) 



where [/(x ± xo,P) are the QC pulses with power P, which are centered at x = T^o? and r]i^2 are the kicks applied 

to them. Simulations of Eq. (|2]) with initial conditions (|22]) are displayed in Fig. [8l 

The simulations demonstrate the repulsive interaction between the colliding QC. In Figs. 8(a) and |8(b)[ they 
bounce back without overlapping, because the collision velocity is small (nevertheless, one of the solitons suffers a 
conspicuous perturbation in Fig. 8(a)[ ). In Figs. |8(c) and |8(d)] the larger velocities give rise to the overlap between 
the colliding QCs, which results in the strong inelasticity of the collision in the case displayed in Fig. |8(d)[ 



IV. CONCLUSIONS 



The objective of this work was to consider ID solitons in the INPhC (inverted nonlinear photonic crystal) with 
competing LL and NL (linear and nonlinear lattices). Unlike recently studied models with the cubic nonlinearity, we 
here consider the saturable self-focusing. The ID crystal of this type was recently fabricated, using the SU-8 polymer 
material periodically doped with Rhodamine B, which lends the medium the saturable nonlinearity. Combining 
numerical methods and analytical approximations, we have demonstrated that broad solitons are sharply localized in 
this setting, thus taking the shape of the QCs (quasi-compactons). With the increase of the total power, P, the QC 
remains a stable object, which switches its position between the linear and nonlinear layers, which form the INPhC. 
The width of the soliton increases with P, which is a manifestation of the saturable nonlinearity. The large width of 
the QC makes it a mobile object, unlike solitons in the INPhC with the cubic nonlinearity. The threshold value of the 
transverse kick (phase tilt), which is necessary to set the QC in motion, was found as a function of P. The threshold 
value gradually decays with the increase of P. Collisions between two moving QCs were are also studied, by means 
of direct simulations. 
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(c) (d) 

FIG. 7: (Color online) (a) Direct simulations of the evolution of the soliton with total power P = 80 in the model with the 
cubic nonlinearity [Eq. (|2T]) ]. initiated by the kick with strengths 77 = 0.15(7r/(i) (a), 0.25(7r/(i) (b), and 77 = 0.3b{7r/d) (c). 
(d) The top plot: ?7ic is the border between stable solitons and fuzzy beams; 772c is the border between the fuzzy beams and 
completely destroyed ones. The bottom plot shows the average width of the quiescent soliton versus its total power P. 



The results reported in this work may be applied to the design of all-optical data-processing schemes. In particular, 
the power-controlled switch of the spatial soliton between adjacent layers, as well as the high mobility of these solitons, 
may be quite relevant properties, in this context. 

This work may be extended in other directions. It particular, 2D waveguide arrays, based on the RhB-SU-8 
mixtures, can be fabricated [U, [22|, which suggests to consider spatial solitons in two-dimensional INPhC with the 
saturable nonlinearity. 

Y.L. appreciates useful discussions with Dr. M. Feng. This work was supported by the Chinese agencies NKBRSF 
(grant No. G2010CB923204), NSFC(grant No. 11104083) and CNNSF(grant No. 10934011). 
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